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Abstract 

Background: Changes in children's contact patterns between termtime and school holidays affect the transmission of 
several respiratory-spread infections. Transmission of varicella zoster virus (VZV), the causative agent of chickenpox, has also 
been linked to the school calendar in several settings, but temporal changes in the proportion of young children attending 
childcare centres may have influenced this relationship. 

Methods:\Ne used two modelling methods (a simple difference equations model and a Time series Susceptible Infectious 
Recovered (TSIR) model) to estimate fortnightly values of a contact parameter (the per capita rate of effective contact 
between two specific individuals), using GP consultation data for chickenpox in England and Wales from 1967-2008. 

Results:lhe estimated contact parameters were 22-31% lower during the summer holiday than during termtime. The 
relationship between the contact parameter and the school calendar did not change markedly over the years analysed. 

Conclusions: In England and Wales, reductions in contact between children during the school summer holiday lead to a 
reduction in the transmission of VZV. These estimates are relevant for predicting how closing schools and nurseries may 
affect an outbreak of an emerging respiratory-spread pathogen. 
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Introduction 

Chickenpox is caused by varicella zoster virus (VZV), which is 
spread primarily by the respiratory route. In temperate regions, 
almost all individuals become infected with VZV during their 
lifetime, usually during childhood [1]. Several modelling studies 
have found that, similar to measles [2-6], VZV transmission is 
reduced during school holidays compared to termtime [5-7] . Two 
of these studies used historical data, from early 20 th Century 
Denmark [5] or from 1928-1973 in the USA [6]. A third study 
used a more recent, but relatively short, time series from France 
(1991-96) [7]. The importance of schools in the transmission of 
VZV is further supported by data from Massachusetts, USA 
(1952-61), which showed that 62% of cases were in children aged 
5—9 years and 19% in children aged 6 years, i.e. having only 
recently started school [8]. 

Since the mid 1 980s, the age distribution of GP (primary care) 
consultations for chickenpox in England and Wales has changed, 
with an increasing proportion of cases being aged <5 years and a 
decreasing proportion being aged 5-14 years [9]. This changing 
age distribution may be due to changes in contact patterns 
resulting from an increasing proportion of children attending 
preschool childcare centres [9-1 1]. This raises the question of how 
much increases in contact between pre-school children have 



reduced the influence of the school calendar on transmission of 
childhood infections. A similar change in the age distribution of 
reported chickenpox cases, attributed to increased urbanisation, 
has been reported for Massachusetts during 1942-51 compared 
with 1952-61, with an increasing proportion of cases being seen in 
children aged under 10 years [8]. 

In this paper, we apply two modelling approaches to long-term 
data (from 1967-2008) from England and Wales to estimate how 
the rate of VZV transmission changes during the year and whether 
its relationship to school holidays has changed over time. We 
calculate fortnighdy values of a contact parameter, defined as the 
per capita rate of effective contact (i.e. contact sufficient to allow 
transmission) between two specific individuals [12]. The estimates 
indicate how school closure during an outbreak of a respiratory- 
spread pathogen might affect contact patterns and transmission of 
infection. 

Methods 

Datasets 

Incidence data. The Royal College of General Practitioners 
(RCGP) Research and Surveillance Centre runs a sentinel primary 
care (GP) surveillance scheme in England and Wales [13]. 
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Currently, data are collected electronically from ~ 1 00 general 
practices in England and Wales (101 practices contributed data in 
2009, covering a population of over 900 000 [14]). The RCGP 
provided weekly overall and age-specific GP consultation rates 
(1967-2008) for chickenpox (Figure SI in File SI). 

Birth rates. Fortnightly birth rates per 100,000, for 1967— 
2008, were estimated by dividing the annual number of births 
(available from the UK Office for National Statistics [15]) by 26 
and using mid-year population estimates for England and W ales as 
the denominator [16—18]. 

School holiday dates. In England and Wales, school holiday 
dates are set locally by 326 Local Authority Districts (LADs). 
School holiday dates were identified from the websites of three 
randomly selected LADs from each of nine geographical regions 
[19-21]. Term dates were also available for the Inner London 
Education Authority (ILEA) for 1952/53 to 1979/80 (ILEA, 
unpublished data), and from the websites of the 13 LADs formed 
from the abolition of ILEA in 1990 [22]. On the basis of these 
data, fortnights 1, 7-8 (Easter, in April), 15-18 (late July to early 
September) and 26 (Christmas) were treated as holidays. 
Sensitivity analyses treated fortnights 1, 8, 16—18 and 26 as 
holidays. 

Modelling methods: overview 

We used two approaches to explore how the rate of effective 
contact differs between termtime and holidays. Firstly, we used a 
simple mass action model, based on one used by Fine and 
Clarkson [2,23], to calculate the contact parameter in each 
fortnight for individual years. In the second approach, we used 
"Time Series Susceptible Infectious Recovered" (TSIR) modelling 
[3], in which the contact parameter can vary seasonally, with a 
pattern assumed to be the same each year. Both methods are 
based on the following equations, which describe the number of 
susceptible (S t ) and infectious (/,) individuals at each time point, t 



S, + B, 



i =P,I,S, 



(1) 



(2) 



/ ;H 1[ represents the number of infections that are reported at time 
t+l; dividing this by the reporting fraction, p (the proportion of 
infections which are reported), provides an estimate of the total 
number of infectious individuals present in the population. B, 
represents the number of births which occurred at time t. The 
equations assume that individuals contact each other randomly. 
The time step used is the serial interval (the time between 
successive cases in a chain of transmission [24]), assumed to be two 
weeks [25]. 

Estimating the reporting fraction 

We used an approach from TSIR modelling to estimate the 

reporting fraction (the percentage of infections that are reported to 

the surveillance system) for both the simple mass action model and 

the TSIR model. 

To estimate the reporting fraction (p), we first plotted the 

i 

cumulative number of infections reported by a given time ( 1^ ) 

k=\ 

against the corresponding cumulative number of births 
t 

( Bk-d), where d represents the duration of maternally derived 



immunity [3], assumed to be six months [26]. We then fitted a 
straight line through this plot and the reciprocal of the gradient of 
this line was taken to equal the reporting fraction. This method for 
estimating p follows from the observation [3] that, if all individuals 
become infected in their lifetime and infection is never fatal, 
Equation 1 can be rewritten as the equation of a straight line 
relating the cumulative number of births to the cumulative 
number of reported infections, with gradient l/p: 



Do 



(3) 



D t represents the residuals from the fitted line, corresponding to 
the difference between the mean number susceptible (S) and the 
number susceptible at time / (i.e. D t = S,— S) [3]. 

As the age distribution of chickenpox consultations changed 
over time [9,11,27] (Figure 1) and symptom severity varies with 
age, we allowed the reporting fraction to differ between three 
periods defined by the age distribution of the consultations, i.e. we 
fitted separate lines to the plot of the cumulative number of cases 
against the cumulative number of births for the following periods: 
1967-76 (most consultations were for 5-14 year-olds), 1977-97 
(the proportion of consultations in 5-14 year-olds declined and 
that in 0-4 year-olds increased) and 1998-2008 (most consulta- 
tions were for 0-4 year-olds). 

The simple mass action model 

We used a modification of the method used by Fine and 
Clarkson [2] [23] to estimate fortnightly values of the contact 
parameter (fi t ) in each year using Equations 1 and 2. Equation 2 
can be rearranged as follows in order to estimate fortnightly values 
of the contact parameter: 



it+i 
i,s, 



(4) 



The values of and were taken from the RCGP data, 
corrected for under-reporting using the reporting fractions 
estimated from the TSIR model. The number of susceptible 
individuals at the start of 1967 was assumed to be either that 
estimated from the TSIR model (described below), the mean of the 
estimated values from the TSIR model when three time periods 
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Figure 1. Proportion of chickenpox consultations in different 
age groups, 1967-2008. 

doi:10.1371/journal.pone.0099762.g001 
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were considered separately (see below for further details), or 20% 
of the population (consistent with the observed proportion of the 
population that was found to be seronegative, according to 
serological data from England and Wales collected in 1996 [1]). 

TSIR model 

Overview of TSIR modeling. The TSIR model [3,5] is an 
extension of the simple mass action model in which Equation 1 is 
rearranged in order to estimate the reporting fraction as above, 
and a linear expression derived from Equation 2 is used to estimate 
fortnightly values of the contact parameter. Unlike the method 
used here with the simple mass action model, the TSIR method 
assumes that any seasonal patterns in the contact parameter are 
the same in all years analysed. Analyses of data using TSIR models 
often try to replicate the predictions that might be obtained in 
models which assume non-homogeneous mixing, by re-writing 
Equation 2 using an additional parameter, a, as: 

It+i=P,ItS, (2) 

a is less than or equal to 1 , with a = 1 corresponding to random 
mixing [3]. Although the interpretation of this parameter is not 
straightforward, its inclusion can help to correct inaccuracies 
resulting from modelling a continuous process using discrete time 
steps [28]. 

Estimating the contact parameters. The contact param- 
eter /?, for a given week t over different time periods (see above) 
and the parameter a were estimated by linear regression using 
Equation 5. This uses the estimated true numbers of chickenpox 
cases per 100,000 in fortnight tior I, (the reported number of cases 
per 100,000 corrected for under-reporting): 

ln7, + i = lnft +<x\nl t + ln(S+ D,)+ e, (5) 

This equation is obtained by substituting S, = S+D t into Equation 2 
and taking logs of the resulting equation; e, t is an error term. 
Marginal profile likelihoods were used to estimate S(see below) [5] . 

This regression was fitted separately for the periods 1967-76, 
1977—97 and 1998-2008. In each regression, variations in ft, 
during the course of the year were assumed to be the same in each 
year (e.g. the contact parameter in fortnight 1 of 1967, 1968, 
1969,. ..1976 was assumed to be identical). The model was fitted 
repeatedly, each time using a different value of S and estimating 
the associated values of/?, and a. Goodness of fit was assessed using 
the log-likelihood deviance between the natural log of the 
estimated number of infections (the reported number corrected 
for under-reporting) and that predicted by Equation 5 as the 
goodness of fit criterion (see the Supplementary Information for 
the definition of the log-likelihood deviance). Values of S used in 
this fitting ranged from the minimum value which would allow S + 
D t to be positive to 100,000 (Equation 5). The best-fitting values of 
ln(/? ( ) and hence the contact parameter in fortnight t (/?,) were taken 
from the model with the value of S which resulted in the smallest 
log-likelihood deviance. 

In sensitivity analyses we repeated the analyses after fixing the 
parameter a to equal 1 and used the Akaike's Information 
Criterion (AIC) [29] to compare the resulting fit against that 
obtained without restricting the size of a. In addition, we 
explored the effect of assuming that the contact parameter 
followed the same cycle during the course of the year from 1967 
to 2008. 

Analyses of the consistency between model estimates and 
the observed data. We assessed the consistency between 



predictions based on the estimated values of fi t and the observed 
data by comparing estimates of the number of cases each week 
(taken as the estimated reporting fraction multiplied by the 
exponents of the fitted values from the regression for each fortnight 
(Equation 5) (i.e. ty) and those reported in the RCGP data using 
the correlation coefficient, r. Also, Equations 1 and 2 were 
evaluated using the estimates of f$ t and a to calculate the number 
of infectious cases each fortnight (per 100,000), which were 
compared against the observed data. The assumed number of 
infectious and susceptible individuals in 1967 were based on the 
reported chickenpox rate from the RCGP data (divided by the 
reporting fraction) and the estimated number susceptible. 

Relationship between the contact parameter and school 
holidays 

For the estimates from the simple mass action and TSIR 
models, the percentage difference between /$, during termtime and 
school holidays (all holidays or the summer holiday) was calculated 
using the expression 100 x ((/? term _-/? holiday )//? term _), where /? is the 
mean value of the contact parameter during the given time period. 
95% CIs for the percentage difference were calculated using the 
bias-corrected and accelerated bootstrapping method [30] . 

All analyses were carried out in Stata, version 12. 

Results 

From 1967-2008, chickenpox consultation rates were highest in 
0-4 and 5-14 year-olds (Figure SI in File SI). Consultation rates 
were highest in spring and lowest during the summer (typically 
reaching a minimum between weeks 37 and 39). 

Estimates of the contact parameters using the simple 
mass action model 

The contact parameter estimated using the simple mass 
action model was in general lowest during the summer holiday, 
but seasonal patterns differed between years (Figure 2). Con- 
sidering all school holidays across all years, the mean value of 
the contact parameter was 11% (95% CI —6-21%) lower 
during holidays than during termtime. Estimates for individual 
years ranged from a 30-35% reduction to a 7-8% increase. 
Most of the point estimates suggested that the contact parameter 
was lower during holidays than termtime, although the 95% CI 
included zero in most years (Figure 3A). Considering just the 
summer holidays, the mean estimate of the contact parameter 
during the holiday was 22% lower than that during termtime 
(95% CI 13-31%, Figure 3B). All of the point estimates for 
individual years indicated a reduction in the contact parameter 
during summer holidays, although this was sometimes small 
(ranging from 0.3% to 45%) and the CIs sometimes included 
zero (Figure 3B). 

Estimates from the simple mass action model of the percentage 
difference between the contact parameter between termtime and 
school holiday were insensitive to the assumed number of 
susceptible individuals in 1967 and the dates of school holidays 
(Figures SI 0-1 2 in File SI). 

Analyses using TSIR modelling 

The estimated reporting fraction was highest during 1977-98 
(42.8%, 95% CI 42.7-42.9%), compared to 29.6% (95% CI 29.3- 
29.9%) and 31.6% (95% CI 31.4-31.8%) in 1967-76 and 1998- 
2008, respectively (Figure S2 in File SI). 

Analyses of the RCGP data for the periods 1967-76, 1977-97 
and 1998-2008, using Equation 5, suggested that on average 9% 
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Figure 2. Fortnightly estimates of the contact parameter for chickenpox based on the simple mass action model. The proportion 
susceptible at the start of the time series was assumed to be 13%; shading shows approximate timing of school holidays. 
doi:1 0.1 371 /journal.pone.0099762.g002 



(95% CI 2-100%), 11% (95% CI 2-100%) and 18% (95% CI 1- 
100)%, respectively, of the population was susceptible to VZV, 
although the CIs were very wide. The corresponding estimates of a 
were 0.728 (95% CI 0.643-0.814), 0.682 (95% CI 0.619-0.744) 
and 0.704 (95% CI 0.617-0.791), respectively. 

Estimates of the contact parameter decreased over time, ranging 
from 2.25xl0~ 4 to 4.71 xl0~ 4 per fortnight in 1967-76, 
2.00xl0" 4 to 4.00x" 4 per fortnight in 1977-97, and 
1.15xl0" 4 to 2.30xl0~ 4 per fortnight in 1998-2008 (Figure 4). 
The mean value of the contact parameter estimated for the 
summer holiday was 31% (95% CI 26-38%), 27% (95% CI 11- 
36%) and 27% (95% CI 8-35%) lower than that during termtime 
in 1967-76, 1977-97 and 1998-2008, respectively (Figures 4 and 
5). All holidays considered together had no clear effect on the 
contact parameter. These estimates were fairly insensitive to the 
assumptions regarding the timing of holidays (Figure 5 and Figure 
S3 in File SI) and were similar to those obtained when a was 



constrained to equal 1 or when the TSIR model was applied to all 
the data from 1967-2008 (Figures S3-S9 in File SI). 

The fitted values of the number of chickenpox cases from the 
regression (scaled down for under-reporting) compared well with 
the observed RCGP data: the correlation coefficient was 0.929 
(95% CI 0.920-0.937, Figure 6A). The output obtained after 
substituting the estimated contact parameters and a into the 
difference equations (Equations 1 and 2) also correlated reasonably 
weU with the data (r= 0.631, 95% CI 0.594-0.666), although the 
correspondence was relatively poor during 1998-2003 (Figure 6B). 

Discussion 

We estimate that contact between children sufficient for 
transmission of chickenpox is 22-31% lower during school 
summer holidays than during termtime. The effect of the summer 
holiday does not appear to have changed substantially over time, 
despite increases in the proportion of preschool children attending 
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Figure 3. Percentage reduction in the contact parameter for chickenpox during school holidays, by year, 1967-2008. Estimates are 
based on the simple mass action model. Error bars show 95% confidence intervals. 
doi:1 0.1 371 /journal.pone.0099762.g003 



childcare [31]. The results based on the TSIR and simple mass 
action models were similar. 

We used a long, contemporary time series to analyse the effects 
of school holidays on contact patterns, taking into account 
temporal changes in the reporting fraction. However, the 
comparisons of the contact parameter between termtime and 
holidays have some limitations. They do not account for other 
seasonal factors, such as weather, which might affect transmission 
(particularly temperature: several studies have found the highest 
incidence of chickenpox to occur in the coolest months of the year) 
[32]. Some fortnights might have been misclassified as holidays, 
due to slight geographic and temporal variations in holiday dates. 
The two-week time step used (consistent with the serial interval of 



chickenpox [25]) also introduces some error as some fortnights 
may comprise one week each of holiday and termtime. Finally, the 
analyses do not account for the possibility that children might be 
more likely to see a GP during termtime than during school 
holidays (e.g. if parents require advice regarding how long their 
child should be excluded from school). This would reduce the 
reporting fraction during holidays compared to termtime, acting to 
overestimate the effects of holidays on the contact parameter. 

We did not formally assess the effects of the Christmas and 
Easter holidays on contact patterns, for several reasons. Firstly, 
consultation behaviour over the Christmas and New Year periods 
will differ from that during the rest of the year. Secondly, the 
Christmas holiday lasts only two weeks, meaning that any 
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Figure 4. Fortnightly estimates of the contact parameter for chickenpox as estimated by TSIR modelling. A) 1967-76; B) 1977-97; C) 
1 998-2008; with values of a as shown. Error bars show 95% confidence intervals. Shaded rectangles show the approximate timing of school holidays. 
Fortnight 1 is the first two weeks of January; fortnight 26 is the last two weeks of December. 
doi:1 0.1 371 /journal.pone.0099762.g004 



comparison of the ji t estimates for termtime and the Christmas 
holiday from the TSIR model would be based on a single value for 
the holiday. Thirdly, the dates of the Easter holiday vary 
considerably between years, so the classification of fortnights 
around this period as holiday or termtime is uncertain. 

The TSIR and simple mass action models are based on the 
same principles, and produced similar point estimates of the effects 
of school holidays on the transmission of VZV. The TSIR model is 
statistically rigorous, e.g. it easily allows calculation of confidence 
intervals for the weekly estimates of the contact parameter. 
However, it assumes that the annual pattern in the contact 
parameter (and therefore the effect of school holidays) is the same 
in all years. Our implementation of the simple mass action model 
allowed the seasonal pattern to differ between years but does not 
indicate the precision of the estimates of the fortnightly contact 
parameters. The absence of a clear trend in the effects of school 
holidays on the contact parameter as estimated using the simple 
mass action model therefore supports the conclusion from the 
TSIR model that there was no systematic change in the effects of 
school holidays over time. 
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Figure 5. Estimated percentage difference between the contact 
parameter (per fortnight) for chickenpox during termtime and 
holidays. A) all school holidays; B) summer holidays, with different 
values of a. A positive value represents a reduction in the contact 
parameter during holidays. 
doi:1 0.1 371 /journal.pone.0099762.g005 



Our results are consistent with previous studies of chickenpox. 
Monthly estimates of the contact parameter in New York City 
(1931-60) from an SIR model were lowest (~65% of the mean 
value, i.e. a 35% reduction) during June and July [6]. The weekly 
contact parameter estimated from an SEIR model applied to data 
from France (1991—96) was ~60% lower during the summer 
holiday than during a typical week in termtime [7]. The greater 
size of this reduction compared to our estimate may reflect 
differences in childcare practices, or in the length of the time 
period considered. 

Monthly estimates of the basic reproduction number (R () , the 
average number of secondary infectious cases resulting from a 
single infectious person introduced into a completely susceptible 
population) for chickenpox, from a TSIR model applied to data 
from early 20 th Century Copenhagen, also appeared lowest during 
the summer holiday [5]. For example, Ro was ~11— 14 during 
June and July, and ~25 in September. 

Previous studies have suggested that the contact parameter for 
chickenpox between school-aged children has increased over time 
[10], which could result from increasing childcare attendance [31]. 
In contrast, our estimates from the TSIR model suggested that the 
contact parameter decreased over time. However, the absolute 
values of fi t from the TSIR model are difficult to interpret, as they 
depend on the proportion of the population assumed to be 
susceptible at the start of the time period analysed. The point 
estimates of the average proportion susceptible increased over 
time, albeit with extremely wide CIs, which may explain the 
apparent decrease in the contact parameter. The results from 
the simple mass action model did not suggest a reduction in the 
contact parameter over time. 

The average percentage of the population that is susceptible to 
chickenpox was estimated as 9%, 11% and 20% during 1967—76, 
1977-97 and 1998-2008 respectively, but with considerable 
uncertainty. These point estimates are lower than published 
values of the proportion seronegative in the UK. Of 2091 
specimens (taken from 1—20 year olds in England and Wales in 
1996), 24% were seronegative [1]. Samples from 1-39 year olds 
from Sheffield taken between 1970 and 1992 varied in the 
proportion seronegative, from 42% (in 1974) to 22% (in 1992) 
[33]. However, the representativeness of these remainder blood 
samples from other laboratory tests is unclear, and the absence of 
age structure in the models makes it difficult to directly compare 
the estimates of the proportion susceptible with empirical data. 

The seasonal pattern in the contact parameter was independent 
of the assumed proportion susceptible at the start of each time 
period, and remained consistent over time despite changes in 
preschool childcare attendance. This may be because preschool 
centres often (although not always) close over summer [34], which 
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Figure 6. Assessment of model fit for chickenpox data, treating 
1967-76, 1977-97 and 1998-2008 separately. A) Relationship 
between the RCGP chickenpox consultation rates and the fitted values 
from the regression (scaled down for under-reporting); B) RCGP data 
and the values predicted by the difference equations using the 
estimated contact parameters. 
doi:1 0.1 371 /journal.pone.0099762.g006 

might augment the seasonal pattern compared to the situation 
where childcare attendance is limited. Increasing use of group 
childcare for school-aged children during holidays might reduce 
the influence of holidays on the contact parameter, so that overall 
the seasonal pattern remained relatively unchanged. 

Our estimates of the reporting fraction (approximately 30%, 
43% and 32% for 1967-76, 1977-97 and 1998-2008 respectively) 
are comparable to previous estimates from analyses of RCGP 
chickenpox data [35-37]. Previous analyses of the RCGP 
chickenpox data estimated the reporting fraction to be 36% by 
comparing consultation data from 1967-85 to serological data 
[35] . Other studies have also suggested a reporting fraction of 38% 
[36]. The age-specific true numbers of chickenpox cases occurring 
between 1991 and 2000 have been estimated by calculating the 
age-specific force of infection based on serological data [37]; 
comparing these estimated numbers of infections to the numbers 
of consultations in the RCGP data (reported in the same paper) 
produces an overall reporting fraction of 4 1 % (43 % for 0—4 year 
olds, 33% for 5-14 year olds, and 52% for 15-44 year olds). We 
were unable to estimate age-specific reporting fractions, but future 
work could extend Equation 3 to an appropriate multiple 
regression in order to estimate several reporting fractions. 
Temporal changes in the reporting fraction may reflect age- 
dependence in the propensity to consult a GP: the estimate 
reporting fraction and the proportion of consultations in 15-44 
year-olds were both highest during 1977—97, and VZV infection is 
more severe in older age groups than in children [9,37]. Other 
authors have fitted more complex functions, such as splines, to 
estimate the reporting fraction for other childhood infections, 



rather than the linear regressions which we used (Equation 3) [38]. 
Although this would have allowed us to analyse the full period 
1967-2008 whilst still accounting for changes in the reporting 
fraction over time, we chose to use the linear regressions for 
simplicity and because the data for each time period were well 
approximated by a straight line (Figure S2 in File SI). 

An alternative method of estimating the reporting fraction 
involves comparing age-specific serological data with the cumu- 
lative proportion of a birth cohort reported to have had the 
infection by a given age, as done for measles [2]. This would less 
easily allow the investigation of changes in the reporting fraction 
over time, as repeat serological surveys would be needed. Age- 
specific seroprevalence data based on sera taken between 1966 
and 1992 in Sheffield, UK, have been published [33] but the 
relatively small sample size in some age groups and the unclear 
generalisability to the wider UK population may limit the 
reliability with which these data could be used to assess temporal 
changes in the reporting fraction. 

Contact surveys (in which participants report the number of 
individuals with whom they make contact) [39] provide an 
alternative method of measuring the effects of school holidays on 
contact patterns. For example, children reported 19% [40] or 50% 
[41] fewer daily contacts during school holidays than termtime in 
Belgium and England, respectively. Results from modelling 
studies, such as ours, provide additional evidence to complement 
these estimates, and refer specifically to effective contact (the 
nature of which may be difficult to define in contact surveys). 

Transmission of other viral infections spread by the respiratory 
route, including measles [2] [5,6], mumps [5,6] and influenza 
[42], appears to be reduced during school holidays to an extent 
similar to that which we estimated for chickenpox [2,6,42]. For 
example, the contact parameter for measles in England and W ales 
(pre- vaccination) was approximately 27% lower during holidays 
than termtime [2], despite the fact that measles is considered to be 
more transmissible than VZV [43]. Influenza surveillance data 
from France (1984-2006) suggest a 20-29% reduction in 
transmission between children during school holidays [42]. 
Transmission is likely to be affected by school closures only if 
they last longer than the serial interval of the given infection [9]. 

School closures, with or without concurrent closure of nurseries, 
may be considered as an outbreak control measure [44,45], 
particularly if pharmaceutical interventions are unavailable. Our 
estimates of the effect of school holidays on transmission of 
chickenpox are indicative of how such closures might influence 
effective contact patterns for a respiratory-spread infection. 

Supporting Information 

File SI Supporting information, figures and tables. This 
file contains: Calculation of the deviance, Alternative TSIR 
modelling analyses, Additional analyses using the simple mass 
action model, and Figure SI -Figure SI 2. Figure SI, Weekly 
overall and age-specific GP consultation rates for chickenpox, 
England and Wales 1967-2008 (RCGP). Figure S2, Plot of the 
cumulative birth rate per 100,000 against the cumulative reported 
rate of chickenpox per 100,000, 1967-2008. Figure S3, Estimated 
percentage difference between the contact parameter (per 
fortnight) for chickenpox during termtime and that during A) all 
school holidays (fortnights 1, 8, 16-18 and 26) and B) summer 
holidays (fortnights 16-18), as estimated from the TSIR model. 
Figure S4, Fortnightly estimates of the contact parameter for 
chickenpox for 1967-76 (column 1), 1977-97 (column 2) and 
1998-2008 (column 3), assuming a = 1. Figure S5, Assessment of 
model fit for chickenpox data, assuming that a = 1 . A) Relationship 
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between the RCGP chickenpox consultation rates and the fitted 
values from the regression for the periods 1967-76, 1977-1997 
and 1998-2008 separately; B) Time series of RCGP data and 
predictions from difference equations for the periods 1967—76, 
1977-97 and 1998-2008 separately; C) and D) as A) and B) but 
based on the full period 1967-2008. Figure S6, Fortnightly 
estimates of the contact parameter for chickenpox (1967-2008), 
for a = 0.717. Figure S7, Assessment of model fit for chickenpox 
data, using the full period 1 967-2008 and estimating a as 0. 7 1 7 . A) 
Relationship between the RCGP chickenpox consultation rates 
and the fitted values from the regression; B) Relationship between 
the RCGP data and the values predicted by the difference 
equations using the estimated contact parameters. Figure S8, 
Estimated percentage reductions in the contact parameter (per 
fortnight) for chickenpox during school holidays, based on the full 
period 1967-2008. Figure S9, Fortnightly estimates of the contact 
parameter for chickenpox, based on the full time period 1967- 
2008 and assuming a = 1. Figure S10, Fortnighdy estimates of the 
contact parameter for chickenpox based on the simple mass action 
model and RCGP consultation data, assuming that 20% of the 
population was susceptible to infection at the beginning of the time 
series. Figure SI 1, Fortnightly estimates of the contact parameter 



for chickenpox based on the simple mass action model and RCGP 
consultation data and assuming the proportion susceptible at the 
start of the time series to be 12% (the mean of the mean 
percentage susceptible during 1967-76, 1977-97 and 1998-2008); 
13%. Figure S12, Percentage reduction in the contact parameter 
for chickenpox during school holidays, by year, 1967-2008 (simple 
mass action model), based on alternative assumptions about the 
dates of school holidays. 
(DOCX) 
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